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Models where the dark energy is a scalar field with a non-standard Dirac-Born-Infeld (DBI) kinetic 
term are investigated. Scaling solutions are studied and proven to be attractors. The corresponding 
shape of the brane tension and of the potential is also determined and found to be, as in the 
standard case, either exponentials or power-law of the DBI field. In these scenarios, in contrast to 
the standard situation, the vacuum expectation value of the field at small redshifts can be small 
in comparison to the Planck mass which could be an advantage from the model building point of 
view. This situation arises when the present-day value of the Lorentz factor is large, this property 
being per se interesting. Serious shortcomings are also present such as the fact that, for simple 
potentials, the equation of state appears to be too far from the observational favored value — 1. 
Another problem is that, although simple stringy- inspired models precisely lead to the power-law 
shape that has been shown to possess a tracking behavior, the power index turns out to have the 
wrong sign. Possible solutions to these issues are discussed. 

PACS numbers: 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 



Since the discovery that the expansion of the Universe 
is presently accelerated 0, H 1 H i, i, 0, S, ©, E3, 
various suggestions have been made in order to explain 
this observational fact. Amongst them is the hypothe- 
sis of dark energy, a fluid with a negative pressure rep- 
resenting about 70% of the total energy density in the 
Universe. The question of the physical nature of the 
dark energy has, of course, been widely discussed. The 
most natural candidate, still perfectly compatible with 
all the data available, is the cosmological constant. How- 
ever, the difficulty to reconcile the value of A deduced 
from the observations with the value calculated theo- 
retically [ll[ (maybe too naively?) has prompted the 
study of alternatives. Clearly, a simple scalar field, a 
"quintessence" field, is a natural candidate for such an 
alternative [3 El El El 

EH- Amongst all the possibil- 
ities, scalar fields with inverse power-law potentials have 
attracted lot of interest because, in this case, there is 
a solution of the equations of motion that is an attrac- 
tor [l2j • This means that the present-day behavior of the 
Universe is insensitive to the initial conditions. Usually, 
the attractor solution is a scaling solution, i.e. a solution 
for which the energy density scales as a power of the scale 
factor p7| . 

If the above mentioned route is correct, then another 
interesting issue is whether a candidate for quintessence 
in high energy physics can be identified. Clearly, this can- 
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not be done without going beyond the standard model 
of particle physics. In particular, it would be very in- 
teresting to achieve this goal in string theory since it is 
presently our best candidate as a unified theory [l8|. 

Recently, there have been many works aiming at con- 
necting string theory with inflation which is also a phase 
of accelerated expansion (but taking place in the very 
early Universe at a much higher energy scale). For this 
purpose, new ideas in string theory based on the concept 
of branes have revealed themselves especially fruitful. In 
particular, scenarios where the inflaton is interpreted as 
the distance between two branes moving in the extra di- 
mensions along a warped throat have given rise to many 
interesting studies [H, [H HH [Hj]. In this article, we 
want to investigate whether the same kind of ideas can 
lead to sensible dark energy scenarios. 

At the technical level, scenarios of the type mentioned 
above lead to scalar field models where the kinetic term 
is non-canonical. More precisely, the kinetic term has 
a Dirac-Born-Infeld (DBI) form. Physically, this origi- 
nates from the fact that the action of the system is pro- 
portional to the volume traced out by the brane during 
its motion. This volume is given by the square-root of 
the induced metric which automatically leads to a DBI 
kinetic term. Therefore, as a first step toward a sce- 
nario of "DBI-essence" , it is first necessary to understand 
whether scaling and attractor solutions are still present 
when the scalar field has a DBI kinetic term. This ques- 
tion constitutes the main target of the present article. 

This paper is organized as follows. In Sec. [Til we briefly 
review the scaling properties of a quintessence field with 
a standard kinetic term. Then, in Sec. IIII1 we reconsider 
this question but with a DBI kinetic term. In particular, 
we compare the DBI results with the standard ones. In 
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Sec. IIV1 we study the behavior of DBI-essence at small 
redshifts. Since, in this case, the scalar field is no longer a 
test field, this requires numerical computations. Finally, 
in Sec. El we present our conclusions and discuss the 
open issues that should be studied in the future. 



II. SCALING SOLUTIONS WITH A STANDARD 
KINETIC TERM 

We consider a spatially flat Friedmann-Lemaitre- 
Robertson- Walker (FLRW) Universe containing a perfect 
fluid and a scalar field <f>. Assuming that the scalar field 
and the perfect fluid are separately conserved, the equa- 
tions of motion are given by 



p + 3H(p 
+ 3H (p 4 



H 2 
hp) 



3 
0. 

0. 



(p + P<t) 



(1) 
(2) 
(3) 



where H = a /a is the Hubble parameter and p and p are 
respectively the energy density and the pressure of the 
perfect fluid. A dot denotes a derivative with respect to 
cosmic time and the quantity k is defined by k = 8n/m 2 r 
In the following, we assume that the perfect fluid has a 
constant equation of state parameter oj = p/p, the two 
cases of main interest being u> = 1/3 for the radiation- 
dominated era and to — for the matter-dominated 
era. In this case, the conservation equation @ can 
be integrated exactly and leads to the familiar behav- 
ior p oc a~ 3 ' 1+w ). Moreover, if we further assume that 
the scalar field is a test field and that the evolution of the 
background geometry is mainly controlled by the perfect 
fluid, then one has a(t) oc t 2 /[ 3 ( 1+LJ )l or, for the Hubble 
parameter, H = 2/ [3(1 + ui)t]. 

Let us first briefly remind what are the scaling solu- 
tions in the simple case where the scalar field has a stan- 
dard kinetic term. In this situation, the energy density 
and the pressure are given by the familiar expressions 



P<t> 



V(<f>), p 



v{4>) 



(4) 



If one inserts these expressions into the conservation 
equation for the scalar field ([3]), then one obtains the 
Klein-Gordon equation 



4> + 3H<j> + V'(ct)) = 0, 



(5) 



where a prime denotes a derivative with respect to <j). 
Then, we seek potentials V(4>) such that the energy den- 
sity of the test scalar field scales as a power-law of the 
scale factor, namely p^ oc a~ 3 ( 1+ " < ^ where = vsl Pi, 
is a constant. It has been established in Ref. [12| that 
scaling solutions exist if the potential has an exponential 
shape, 



v{4>) 



M 4 e 



(6) 



where A is a constant or is of the Ratra-Peebles type 
(i.e. inverse power-law of the field), namely 



V(4>) = Af 4+a </r" 



(7) 



In the first case, the particular solution leading to the 
scaling behavior reads 



(S) 



the constant A and the mass scale M being linked by the 
relation 



\ 2 mH 2 = 



2(1 -ui) 
l+u> 



(9) 



As is well-known, the quintessence equation of state pa- 
rameter is just given by the equation of state of the back- 
ground perfect fluid, to^ — w. The particular solution (jSJ) 
is important because it is an attractor. This means that 
the final (i.e. present day) evolution of the field is in fact 
independent of the initial conditions. At the technical 
level, this can be seen by studying small (linear) per- 
turbations around the attractor. The eigenvalues of the 
perturbations around the critical point can be expressed 
as 



A+ = 



1 

2m 



(m — 6) ± \J (m — 6) 2 + 8m (m — 6) 



(10) 



where we have defined m = 3(1 + uj). Since w < 1, 
one has m < 6 and the eigenvalues are negative and one 
has a stable spiral point. Moreover, for the particular 
solution (El), one has 



d^ = 9 

dcj) 2 2 V 



)H 2 



(11) 



This is an important formula because it implies that 
(f> ~ m P1 today. Indeed, V" ~ V/(/) 2 and H 2 ~ V/m^ 
when the field starts dominating the energy density con- 
tent of the Universe; equating these two quantities leads 
to the above mentioned conclusion. For this reason, a 
sensible model building of quintessence is only possible 
in a supergravity (SUGRA) framework [3, HE]. How- 
ever, a well-known difficulty of the exponential case is 
that the property lo^ = lu implies that the scalar field 
cannot drive an accelerated expansion. This is why the 
inverse power-law case seems to be more interesting. 

In the case of the Ratra-Peebles potential ([7]), there 
also exists an exact particular solution of the Klein- 
Gordon equation that is an attractor. It reads 



2/(a+2) 



(12) 



where the quantity </>o is linked to the mass scale M by 
the formula 



M i+a <f> a - 2 t 2 = 



a(a + 2)V1+^ a + 2 / 



(13) 
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For this particular solution, the equation of state param- 
eter can be expressed as 



aoj - 



u>(j> 



(14) 



As expected, in the limit a — > +00 , one recovers the ex- 
ponential case, lo^ = uj. However, the crucial difference 
with the exponential potential is that one can now have 
lo^ < uj, that is to say, the scalar field energy density can 
now scale more slowly than the background fluid and, 
hence, eventually dominates, causing the Universe to ac- 
celerate. Moreover, the solution (fT2")) is also an attrac- 
tor as revealed by a dynamical system analysis. Indeed, 
the eigenvalues of small perturbations around the critical 
point read 



A± = 



(2n 



6)±J(2n 



6) 2 + 8m (n-6) 



2m 



( 15 ) 

where we have defined n = 3(1 + lu^). Again, one can 
show that there exists a stable spiral point as long as both 
eigenvalues are negative, which is equivalent to 2n — m — 
6 < 0. On this attractor, the evolution of the second 
order derivative of the potential is given by 



d 2 V 



9a 
2~< 



1 



(1 



H 1 



(16) 



and one check that this last equation reproduces the 
corresponding equation in the exponential case when 
a — > +00. Again, this prompts a SUGRA treatment 
of the model building issue since one still has <j> ~ m P1 
now. 



III. SCALING SOLUTIONS WITH A DBI 
KINETIC TERM 

Let us now consider that the dark energy scalar field is 
a Dirac-Born-Infeld (DBI) scalar field. In this case, the 
action of the field can be written as 



S DBI = -J d*xa 3 (t) 



T(0h 1 



V(<f>) - T{4>) 



(17) 

where T(</>) is the tension and V(<p) is the potential. From 
this expression, it is easy to obtain the corresponding 
energy density and pressure of the scalar field. They 
read 



^ = ( 7 -l)T(0) + y(<£>, 



7-1 

7 



T(</>) - V(cj>) , 
(18) 

where the quantity 7 is reminiscent from the usual rela- 
tivistic Lorentz factor and is given by 



1 



7 



1 - 03/ T (0) 



(19) 



The expressions (fTHJ) of the energy density and pressure of 
the DBI field should be compared to their standard coun- 
terpart, see Eqs. (j4|). As usual, if one inserts Eqs. (fT8|) 
in the conservation equation ([3]), one obtains the DBI 
Klein-Gordon equation, namely 



3T'M 
2T(cj>) ' 



T'(0)]=O. 
(20) 

We notice that the equation of motion for cf> is quite com- 
plicated compared to Eq. ([5]) despite the fact that the 
conservation equation has retained its standard form. 

Let us also compare with other works in the litera- 
ture. Let us start with K-essence where the action can 
be written as 



d 4 x. 



-9P(4>,X) 



(21) 



where X = (V(f>) /2. Clearly the action (|17p is a spe- 
cial case of the above action. However, as first dis- 
cussed in Ref . [23| , K-essence usually means that the po- 
tential term vanishes and the negative pressure of the 
scalar field is realized only by considering the kinetic 
term [U, [H, [H, H3] . On the other hand, our model can- 
not realize the negative pressure without the potential 
term, as shown below. Therefore, this class of models 
cannot encompass Eq. (TIT)) . Another model related to 
the present study is the case where the dark energy field 
is a tachyon for which the action is given by 28] 



S = - 



(22) 



where M is a fundamental scale and V(T) is a poten- 
tial which, of course, needs not to be the same function 
as V((f>) in Eq. (jTTJ) . This class of theory is equivalent 
to the case studied here (through a redefinition of the 
field) only when the potential in Eq. (TP?]) vanishes. Let 
us also notice that when V(T) is constant the model is 
in fact equivalent to the Chaplygin gas with the equa- 
tion of state p cx — 1/p [29|, [3(| ■ Therefore, beside the 
fact that the search of scaling solutions has not yet been 
investigated in this type of models, we conclude that the 
class of scenarios under scrutiny in this paper was not 
considered before. 

As a warm up, let us find the scaling solutions in the 
simple case where the potential vanishes. As already 
mentioned before, this means that we now seek tensions 
T(4>) such that the energy density of the test DBI scalar 
field scales as a cx a~ 3 ( 1+UJ,t, \ From Eqs. (fTSJ), it is 
easy to show that 7 is, in this case, constant and given 
by 7 = l/w^. Then, the formula of the energy density, 
/?0 = (7 — 1)T{4>) cx a _3 ' 1+w *^ immediately gives the 
scaling in time of T((f>) which in turn, combined with 
4> 2 /T{4>) — (7 2 — l)/7 2 and the fact that 7 is constant, 
implies that 



<f> cx t 



-(i+^)/(i+w) 



(23) 
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This equation is easily solved. Let us start with lu^ = u>. 
In this case, one has <f> oc lni and, as a consequence, 



T{4>) = M 4 



-\<t> 



(24) 



where A is a constant. Again, this case is very simi- 
lar to the situation where we have a standard kinetic 
term and a exponential potential. As a consequence, 
this model suffers from the standard phenomenological 
problems. Since the scalar fields exactly tracks the back- 
ground matter, one cannot have a large enough contri- 
bution of dark energy density today without spoiling Big 
Bang Nucleosynthesis (BBN). Moreover, the scalar field 
behaves as matter today and, therefore, cannot cause the 
acceleration of the Universe. 

On the other hand, if uj^ ^ lo, then the scalar field is 
just a power-law of the cosmic time which implies that 
T{4>) can be expressed as 



T(4>) = M i+a 4>- a 



(25) 



where M is a mass scale and lo^ is related to a and the 
background equation of state parameter uj through the 
relation 



1 

1 



(26) 



Interestingly enough, as discussed before, this is exactly 
the equation obtained when there is a standard kinetic 
term with a Ratra-Peebles potential, see Eq. (|14p. How- 
ever, in the present case, — 1/7 > 0. This means that 
the solution is physically relevant only if u> > 2/ a which 
excludes the case w = 0, at least for a > 0. 

Since it appears that the previously described situation 
is not satisfactory, we now envisage the case where the 
potential V(<f>) is non vanishing. In order to deal with 
this problem, we assume that 7 is a constant. Without 
this hypothesis, the problem is technically very compli- 
cated but the actual convincing argument in favor of this 
assumption is that the corresponding scaling solutions 
(with 7 constant) are attractors, see below. Then, the 
crucial observation is that Eq. ([23]) is still valid because, 



in its derivation, one has never assumed that V — 0. This 
implies that, as in the case of a vanishing potential, scal- 
ing solutions exist for tensions T(0) given by Eq. (|24|) or 
Eq. (|25[) . Then, the Klein-Gordon equation (|2H)l can be 
used to determine the potential. Straightforward manip- 
ulations lead to 



V(<P) 
T(<f>) 



7 2 -l 

7 



1 



7 



1 + ^0 1 + 7 



(27) 



that is to say the potential is proportional to the tension 
and has also the exponential shape or inverse power-law 
shape. It is interesting to notice that, when the field is 
on tracks in the standard kinetic case, the potential term 
is also proportional to the kinetic term, that is to say the 
ratio of the potential term to the kinetic term K = <j) 2 /2 
is a constant given by V(<j>)/ 'K{4>) — 2/(1 + u^) — 1. 

In the case of an exponential potential, the exact so- 
lution (from now on, we put a subscript "e" to de- 
note the quantities that are evaluated with the exact 
particular solution of the Klein-Gordon equation) reads 
e (i) = 2/Aln(t/*o) with \ 2 M 4 t 2 = 4 7o 2 /(7 c 2 - 1) and the 
equation of state parameter is uj^, = u. In the Ratra- 
Peebles case, one as <j> e (t) = 4>o(t/to) 2 ^ a+2 ' 1 with 



M 



4+a ,-a-2.2 



47c 2 



(a + 2)2( 7 2-l) 



(28) 



and the equation of state has the standard form given by 
Eq. (fT4|) . It is importance to notice that, because we deal 
with a modified Klein-Gordon equation, the expressions 
of \ 2 MHl and M 4+q 0q a ~ 2 t\ are different from the ones 



obtained previously, see Eqs. (|9|) and fjl3|) . Let us also 
remark that these formulas can either be obtained from 
the requirement that the Lorentz factor is constant or by 
brute force calculation using the Klein-Gordon equation. 

Let us now study the behavior of small perturbations 
around the particular solutions. Let us first start with 
the exponential case. For this purpose, we rewrite the 
equation of motion in terms of u(t) defined by u = \{<j> — 
4> c ) and t = e T . If we write p = u', then one obtains the 
system 



dp 

d7 



1 1 



l + w V7 7o 



A 2 M 4 (e- u - 1) 



AM 



7e 2 - 1 ( 1 



1+LO 1 + 7c 



- 1 



7° 7e 



du 

7w 



(29) 
(30) 



where, now, the quantity 7 is no longer a constant and can be written as 



1 - 



7 2 - 1 



e" 1 + r 



-1/2 



(31) 
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As a consequence, we see that the critical point is (p, u) = 
(0,0). Notice that, for the critical point, one checks that 
7 = 7<j. We now consider the behavior of small per- 
turbations (Sp,8u) around the critical point (0,0). It is 
straightforward to establish that 



_d_ f5p 

dr \Su 



1 



2 




(32) 



Then, the eigenvalues of this matrix can be expressed as 



(m — 6) 



±\l(m-6) 2 +8m (m-3 7 ° ^ 



(33) 



This expression should be compared with Eq. (|10p . The 
only difference is the presence of the factor 7 C in the last 
term inside the square root. Otherwise, and this is quite 
remarkable, the expression is the same. Let us also notice 
that the condition that the kinetic energy cannot exceed 
the total energy is equivalent to n < 3(7 + l)/7c- Since 
7o > 1, the condition n < 3(j 2 + l)/7e is stronger than 
the condition n < 3(-f c + l)/7c and, therefore, is not 
automatically satisfied in our case. We conclude that 
there is a stable spiral point if n < 3(7^ + l)/7e- 

Let us now turn to the inverse power-law case. This 
time, the dimensionless function u(t) is defined by it(r) = 
4>/4>c the definition of the time r remaining the same. 
Then, a straightforward calculation leads to the following 
system of equations 



dp 

d7 



du 

dT 



(ha + 2 2 1 



\a + 2 7 2 1 + cj 



3a p 2 

P~ -7, — 



2 u l + wa + 2 V7 2 7< 2 



1 1 



p. 



r 7c 2 - 1 



1 + u 1 + 7 C 



- 1 



- - - ) + aM^ a <t>- 

a-l 



4+a^-a-2 / _ a _l_ n 



-71^°' 



(34) 
(35) 



where, this time, the Lorentz factor can be written as 

-1/2 



7 



1 



M 4 + a c/) 



a-1 



V 



Aup 



(a + 2Y 



It is clear form the above system that the critical point 
is now given by (p, u) = (0, 1). The next step is to study 
the behavior of small perturbations (dp, 1 + du) around 
the critical point. One arrives at 



_d_ (bp 
dr \Su 



I »- 2 _ 2 2a 
= a+2 1+uj a+2 

\ 1 




(37) 



As expected, in the limit a — ► +00, the above matrix ex- 
actly reproduces the matrix obtained in the exponential 
case, see Eq. (|32|). Then, the next step is to determine 
the eigenvalues. The result reads 



A± = 
± 



2n — m — 6 
2m 




(2n — m — 6) + 8m I n 



3^i I (.38) 

7o 



This expression should be compared with Eq. (|15p . As it 
was the case before, the modification introduced by the 
DBI kinetic term is only apparent in the last term inside 



the square root. Therefore, there is a stable spiral point 
if n < 3(7^ + l)/7 2 and 2n — m — 6 < are satisfied. 

Finally, on the attractor, in the exponential case, it is 
easy to establish that the following relation holds 



d 2 V 



9 7o (1 + lo) 



l-(l + w) 



7e 



1 + 7o 



(39) 



This formula is the generalization of Eq. (fTl"]) . Obviously, 
one can also establish the corresponding expression in the 
inverse power-law case. It reads 



d 2 V 



a + 1 



7o (1 + uty) 



1 - (l + c^) 



7e 



l+7c 



H' 2 



and this is equivalent to Eq. (|16p . This has important 
consequences for model building. Indeed, if one repeats 
the discussion after Eq. ([TTj) . then one arrives at the con- 
clusion that 

4>~^ (41) 



because the second term in the bracket in Eq. (140(1 cannot 
exceed unity. Therefore, if 7 C S> 1, then the vacuum 
expectation value of the field is not necessarily large in 
Planck units. This is certainly an important advantage 
of the DBI models over the standard ones as respect to 
model building issues. 

Finally, it is also worth commenting about the shape 
of the tension T(</>). From a stringy point of view, the 
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FIG. 1: Evolution of the DBI energy density for different 
initial conditions in the case where the potential is of the 
Ratra-Peebles type with a = 4. The value of C is chosen to 
be C ~ 3.443 which corresponds to 7 C = 20, see Eq. (g2}. The 
initial velocity </>i n j is always chosen such that y m \ = 5. The 
solid line corresponds to an initial vacuum expectation value 
of 0ini/?ripi ~ 10 -10 , the dotted line to 0mi/m P i ~ 10 -9 and 
the dashed line to 0mi/m P i ~ 10~ 8 . The energy density of 
radiation (dotted-dashed line) and cold dark matter (dotted- 
dotted-dashed line) are also represented. 



FIG. 2: Evolution of the DBI equation of state for different 
initial conditions in the case where the potential is of the 
Ratra-Peebles potential with a = 4. As in Fig. [JJ the value 
of C is chosen to be C ~ 3.443 which corresponds to 7 C = 20 
and the initial velocity </>j n j is such that 7i n i = 5. The solid 
line corresponds to an initial vacuum expectation value of 
(piui/mpi ~ 10 -10 , the dotted line to <f)ini/rn P1 ~ 10~ 9 and the 
dashed line to </>i n i/m-pi ~ 10 -8 . The final (present day) value 
of 7 is 70 ~ 3.96 and (j>o/mp\ ~ 1.37. Finally, the equation of 
state is such that ujq ~ —0.42, wi ~ 6.68 x 10~ 2 . 



inverse of T(<p) represents the warp factor of the throat 
in which the branes are living. A natural choice [2l| is 
T((j>) oc 4> 4 , that is to say a = —4. Therefore, this case 
belongs to the class of tracking models considered here 
which is a non trivial result. Unfortunately, the sign of 
the exponent is not the correct one. Indeed, for a = —4, 
one has oj^ = 2uj + 1 > uj which means that, despite the 
presence of an attractor, the scalar field scales faster than 
the background fluid and, hence, can never dominate the 
matter content of the Universe. 



IV. NUMERICAL CALCULATIONS 

In this section, we investigate the behavior of the DBI 
scalar field at small redshifts, when it starts dominating 
the matter content of the universe. In this situation, 
the assumption that it is a test field breaks down and 
numerical calculations are required. 

We first check that the attractor is observed numer- 
ically. As a representative example, we have chosen to 
investigate the case a = 4. In Fig. [TJ we have repre- 
sented the evolution of the DBI energy density for three 
different initial conditions (more precisely, the initial ve- 
locity is always the same and corresponds to an initial 
value of the Lorentz factor 7^ = 5 but different initial 
vacuum expectation values 4>- ln i are considered) . In order 
to have a DBI energy density today equal to 70% of the 



critical energy density, we have tuned the scale M of the 
brane tension T ((/>), see Eq. (|25p. The mass scale of the 
potential is determined by Eq. (|27[) which implies that 
V((j)) = CM i+a 4>- a where C is defined by 



7 2 - 1 

C = 2s — t. 



l 



1 + 1 + 7c 



7e 



a(l+uj) l+7c 



(42) 



Choosing a value of C is in fact equivalent to choosing 
the value of the Lorentz factor on the attractor, 7 e , dur- 
ing a phase of evolution characterized by the background 
equation of state lu. So, for instance, in Fig. [TJ we have 
chosen j c = 20 and u> = 1/3. This means that the at- 
tractor solution should be such that j = 20 during the 
radiation dominated era. Given Eq. (|4"2")l and a = 4, this 
choice implies that C ~ 3.443. In this case, there is also 
an attractor during the matter dominated era but the 
corresponding value of the Lorentz factor is different. It 
is easy to show that it reads 



„cdm 



a(C-l) 




a 2 (C-iy 



(43) 



In the present case, this gives 7c dm ~ 5.438. 

The attractor behavior is clearly seen in Fig. [TJ For 
initial conditions (f>i a i/m Pl ~ 10 -9 and </>ini/?7ipi ~ 10~ 8 , 
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the attractor is joined during the radiation-dominated 
era while for </>j n i/m P1 ~ 10~ 10 , it is reached during the 
matter-dominated era. 

In Fig. [21 we have represented the evolution of the 
equation of state for the same situation. Again, the at- 
tractor behavior is clearly noticed. We can even check 
numerically that, on the attractor, Eq. (| 14[) is valid. 
Since we consider a model with a = 4, the DBI equa- 
tion of state during the radiation dominated era should 
be uj<p ~ —0.11. Clearly, this is what is obtained in 
Fig. [2] The present day value of the equation of state 
is ujq ~ —0.42 (the derivative of the equation of state 
at vanishing redshift being uj\ ~ 6.68 x 10 -2 ). The 
corresponding value for a scalar field with a standard 
kinetic term and the same Ratra-Peebles potential is 
ujq ~ —0.487. Firstly, and contrary to a naive expec- 
tation, the equation of state is not pushed towards — 1. 
Therefore, it seems that we do not gain anything in com- 
parison with the model with a standard kinetic term. 
Secondly, the value obtained seems to be too large given 
the constraints available on ujq . Even if one should put a 
damper on these constraints since they have not been ob- 
tained for the model under considerations here (usually, 
a simple law of the form u> = u>q + u>xZ are used, which is 
clearly not valid for the model under consideration here, 
and this can cause a "bias problem", see Ref. [3l|), the 
value is so far from —1 that the model is probably in trou- 
ble from the observational point of view. This is clearly 
a very serious problem for the class of models studied 
in the present article. One possibility is to decrease the 
value of a. For instance, a = 0.3 implies us^ ~ —0.9. 
Of course, the corresponding model seems contrived and, 
in addition, in this case, a small value of the equation 
of state would also be obtained with a standard kinetic 
term. Another possibility would be to consider other 
shapes for the tension and the potential. The new shape 
of T(ip) and V(<j>) should approximatively reduces to the 
inverse power-law shape at large redshifts such that the 
attractor behavior is preserved and should differ from 
it at small redshifts in order to obtain an equation of 
state closer to —1. A typical example is provided by the 
SUGRA potential [2 EH- However, we show below that 
even this solution could not work in the DBI case. 

Let us now study the evolution of the Lorentz factor 
7. It is represented in Fig. [3] for different initial con- 
ditions, similar to the ones considered in the previous 
figures. In particular, the attractor value, valid during 
the radiation dominated era, j c = 20, is clearly seen 
on this plot. The value 7c dm ~ 5.4, valid during the 
matter dominated era, can also be noticed. An interest- 
ing point is that the present value of the Lorentz factor 
is far from 1. For the case under consideration, it is 
70 ~ 3.96. This means that the non-standard kinetic 
term still plays a role even today. As noticed earlier, this 
can have important implications for model building is- 
sues since a large 7 e implies a small vacuum expectation 
value of the field. We have studied this point in more 
details in Fig. 2] where the evolution of the Lorentz fac- 
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FIG. 3: Evolution of the Lorentz factor for different initial 
conditions for the Ratra-Peebles potential with a — 4. The 
value of C is chosen to be C ~ 3.443 which corresponds to 
7c = 20 and the initial velocity </>i n i is chosen such that 7i n j = 
5. The solid line corresponds to an initial vacuum expectation 
value of 4>xnx/m,pi ~ 10~ 10 , the dotted line to 4>im/rn P1 ~ 10~ 9 
and the dashed line to (feu/fipi ~ 1CP 8 . The final value of 
the Lorentz factor is 70 ~ 3.96 



tor is represented for different values of 7 e . We notice 
that the larger j e , the larger the present day value 70 
and the smaller fa. For instance, for 7 C = 1000, one ob- 
tains 70 ~ 176 and 4>a/m rl ~ 0.2 < 1. This is certainly a 
desirable feature since, usually, vacuum expectation val- 
ues of the order of the Planck mass are at the origin of 
many serious problems as, for instance, a coupling with 
the observable sector which violates the constraint on 
the presence of a fifth force and/or on the weak equiv- 
alence principle [H, [H, [34], HH| ■ On the other hand, if 
the vacuum expectation value remains small in compar- 
ison with the Planck mass, then it could be difficult to 
use the SUGRA potential model to push the equation of 
state towards —1. Therefore, we fac e ag ain the "no-go 
theorem" discussed recently in Refs. [33l l34l [3^| : what 
is interesting from the cosmological point (a large vac- 
uum expectation value in order to have loq close to — 1) 
seems to be incompatible with local tests of gravity (a 
large vacuum expectation value usually means a strong 
coupling with ordinary matter). 



V. CONCLUSIONS 

In this section, we recap our main findings and discuss 
further issues that should be investigated. We have stud- 
ied scenarios where the dark energy is a scalar field with a 
DBI kinetic term. We have shown that, if the brane ten- 
sion and the potential possess either an exponential or a 
power law shape, then there exists scaling solutions that 
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FIG. 4: Evolution of the Lorentz factor for different values 
of 7c in the case where the potential is of the Ratra-Peebles 
type with a = 4. The values considered are 7 e = 20 (solid 
line), 7o = 100 (dotted line), 7 C = 500 (dashed line) and 
7o = 1000 (dotted-dashed line). The corresponding values of 
the Lorentz factor today are 70 ~ 3.96, 70 ~ 17.87, 70 ~ 88.56 
and 70 ~ 176.17 respectively. The final vacuum expectation 
values of the DBI field are 4>o/m P1 ~ 1.37, 4>o/mp\ ~ 0.66, 
4>o/m P1 ~ 0.29 and 0o/rripi ~ 0.2 respectively. 



are attractors. Moreover, if the Lorentz factor is large to- 
day, then the vacuum expectation value of the field can 
be small in comparison with the Planck mass. Let us also 
notice that the fact that the scaling solutions obtained 
in this article correspond to brane tensions of the power- 
law form is fairly remarkable since this is precisely what 
happens in simple string-inspired models. Unfortunately, 
one needs a < for which the dark energy density scales 
faster than that of the background. Maybe, the most 
problematic aspects of the scenario is the fact that the 
equation of state today is too far from — 1. 

In order to improve the above described situation, one 
probably needs more complicated string inspired models. 
In particular, one needs shapes of T((f>) and V(cp) that, 
for <fr <C m m / ^/t^T are of the power law form (with a > 0) 
in order to preserve the attractor, and, for <j> m P1 /^/7e, 
deviate from this form in order to push the equation of 
state towards —1. Let us recall at this stage that this is 
exactly what the SUGRA model does, the characteristic 
scale being the Planck mass instead of m Pl /-y/%. 

The fact that the present-day value of 7 G can be large is 



also an interesting feature of the models under scrutiny. 
For example, this implies that the sound velocity squared 
Cg can significantly deviate from 1 in contrast to the stan- 
dard case. Indeed, in the DBI case, the sound velocity 
squared Cg is given by the following expression 



1 



(44) 



2 = dp ( dp 

Cs dx\dxj x 

where p = p(X,<f>), p = p(X,<f>) (and X = j) 2 /2). 
A dark energy component with « 1 implies less 
power on large scales and, hence, could account for 
the low multipoles of the cosmic microwave background 
anisotropics. Moreover, this would also produce peculiar 
features in the matter power spectrum as discussed in 
Refs. [H, H3, H^|. All these properties could be used to 
distinguish the DBI models from the standard ones. One 
more general grounds, it is clear that a complete calcu- 
lation of the dark energy perturbations could bring new 
insights on the model. 

Finally, another interesting issue is that of the coupling 
of dark energy with the rest of the world. As already 
mentioned, this is usually a problem for quintessence 
because a small mass means a force with a very long 
range, see Refs. [3^, H3, HH| . In some scenarios, this also 
implies variation of the constants, as, for instance, the 
fine structure constant. However, in the present con- 
text, the couplings are totally different. For example, 
the coupling with the electromagnetic field is of the form 
-\/det (g M „ + Fftv), where, here, is the induced met- 
ric on the brane. Therefore, one can maybe expect this 
issue to be less problematic than in the standard case. 
More work is clearly needed in order to draw definitive 
conclusions on these matters. 
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